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00 ' Different non-equilibrium situations have recently been considered when study- 

f - **) ing the thermal Casimir-Polder interaction with a body. We show that the 

i Keldysh Green function method provides a very general common framework 

■ for such studies where non-equilibrium of either the atom or the body with 

the environment can be accounted for. We apply the results to the case of 
ground state polar molecules out of equilibrium with their environment, ob- 
serving several striking effects. We consider thermal Casimir-Polder potentials 
in planar configurations, and new results for a molecule in a cylindrical cav- 
ity are reported, showing similar characteristic behaviour as found in planar 
geometry. 



1. Introduction 

Casimir-Polder (CP) or retarded Van der Waals forces 1 is the name given 
to electromagnetic dispersion forces between electrically neutral, but po- 
larisable particles (atoms, molecules) and macroscopic objects. In the 
present paper we discuss CP potentials on particles when the atom-body- 
environment system is not in thermal equilibrium. There exists a rich lit- 
erature on CP forces in thermal equilibrium (cf. the citations in Ref. 2), 
and non-equilibrium systems of an excited atom inside planar structures 
have been studied for a long time. 3,4 The latter works showed that an atom 
excited to an energy higher than thermal energies will have a spatially 
oscillating component in both force and heating rate. 
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Two different types of non-equilibrium have recently been studied. The 
first study concerned atoms next to a substrate whose temperature differs 
from the environment temperature 5 ' 6 where the present radiation was as- 
sumed to be practically unable to excite the atoms. It was found that the 
imbalance between environment and substrate temperature may lead to 
strong force components whose sign depends on which of the two temper- 
atures is greater. This prediction was subsequently confirmed in an exper- 
iment by the group of Cornell. 7 

The second study considered particles in an equilibrium thermal back- 
ground of bodies and environment at uniform temperature, but prepared 
in an arbitrary superposition of internal eigenstates, and thus not neces- 
sarily at equilibrium with this background 8 (cf. the similar results reported 
in Ref. 9; for an extension to non-uniform temperature environments, see 
Ref. 10). It was demonstrated that even ground state particles are sub- 
ject to an oscillating force component in the presence of macroscopic bod- 
ies at nonzero temperature. While utterly unobservable for atoms (which 
are essentially in their ground state when thermalized at room tempera- 
ture) , there could be some hope of observing and even using this effect for 
molecules, 11 ' 12 which have excited states of very low energy. 

A similar investigation was recently made 13 by applying Keldysh theory 
to the system of two atoms prepared in an arbitrary state, in the presence of 
an external (thermal) electromagnetic field. In the first half of this article, 
we discuss the correspondence of the Keldysh formalism with both of the 
above non-equilibrium theories (Sec. 2). In the second half (Sec. 3), we 
apply the results for a particle in a uniform-temperature environment to 
planar and cylindrical cavities. 

2. Correspondence of thermal non-equilibrium theories 

The recent studies by Antezza, Pitaevskii and Stringari 5 (APS) on the one 
hand and Buhmann and Scheel 8 (BS) on the other hand both deal with non- 
equilibrium situations. However, the theory of APS is concerned with effects 
due to a thermal disequilibrium between different parts of the environment. 
The theory of BS describes an in some sense opposite situation where the 
focus is on a disequilibrium between the atom in a non-thermalized state 
and its environment which is itself in thermal equilibrium with uniform 
temperature. 

In this section we demonstrate how the two situations may be bridged 
through a non-equilibrium theory of CP forces based on the Keldysh Green 
function method 13 by showing that the known results for the BS and APS 
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configurations can both be obtained within this formalism. We restrict our 
derivation to the interaction between a ground-state particle and a non- 
magnetic medium embedded in an external (thermal) electromagnetic field, 
neglecting the non-equilibrium dynamics. A more general situation will be 
considered elsewhere. 

The CP potential of a polarisable ground-state particle at position r 
next to a non-magnetic body embedded in an electromagnetic field can be 
calculated with the help of the Keldysh Green function technique. We calcu- 
late the density matrix of the atom using the Keldysh method. 14 Resticting 
our considerations to the potential at initial time, we may neglect terms 
associated with atomic non-equilibrium dynamics. The shift of the atomc 
ground-state due to the interaction with the electromagnetic field can then 
be found in terms of the photon Green functions (for details, see Ref. 13]: 



U(r) =-^J d " ^'Trlrn |a(w) ■ [%,G (1) (r, r, u) - p(r, r, w)] } , (1) 



where ol(uj) is the ground-state polarizability of the atom and G*- 1 ' (r, r, u) 
is the scattering part of the Green tensor of the electromagnetic field, 
G(r, r',u>), which in turn is the unique solution to 



V x V x 



G(r,r» = 5(r-r') 



(2) 



together with the boundary condition at infinity. Finally, p(r, t;r',t') — 
— i(A(r', t')A(r, t)) T , where A is the vector potential of the external elec- 
tromagnetic field. Note that the Keldysh result offers a clear distinction 
between the influence of the zero-point fluctuations of the field, the first 
term in Eq. (1) and the contributions from external fields as contained in 
the second term. 

Let us first compare the Keldysh method with the BS calculations. 8 
If the environment containing the body and the external field is at global 
thermal equilirium with a uniform temperature Tg, the density matrix p 
can be calculated with the help of the fluctuation-dissipation theorem: 



p(r,r',w) = -2# %,T £ )ImG(r,r',w), 



(3) 



where N(u>, Te) = [cxp(huj j RbTe) — 1] 1 denotes the thermal photon num- 
bers. Substituting this into Eq. (1) and discarding the position-independent 
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contribution associated with the bulk Green tensor (G G^), we find: 



+ 



hfj, 



U{ r )=-^p- / duju 2 [2N(w,T E ) + l}Trlm a(u) ■ G (1) (r, r, 

Z7T 



n Jo 



duuj 2 N(uj,T E )Tr flma(u) • ReG (1) (r, r,w) 



(4) 



The first term in Eq. (4) can be cast into an alternative form by writing 
Imz = (z — z*)/(2i), using the identities G*(w) = G(— u>), ot*(u>) = ot(—u>) 
and making the substitution uj n- — u>. The emerging integral over the entire 
real frequency axis can be completed to a closed contour by adding a van- 
ishing integral over an infinite semi-circle in the upper half of the complex 
frequency plane. Evaluating the contour integral via Cauchy's theorem, we 
are left with the contributions from the poles i£ m of [2N(uj,Te) + 1], viz. 



U(r) =^ k B T J2 'dTr [a(^ m ) ■ G« (r, r, i£ m ) 



m=0 



ft/J-0 



+ — / diouj 2 N(uj,T E )Tr Ima(u) • RcG (L > (r, r, w) 
n Jo 1 



(5) 



with £ m — l-Kvnk^T jh. The prime at the Matsubara sum indicates that the 
m = term carries half-weight. For an isotropic particle, the ground-state 
polarizability in the perturbative limit can be given as 



a(tj) — lim -!- 



lOfcl 



(6) 



[wfe = — E )/h are transition frequencies; d fc are electric dipole matrix 
elements; I is the unit tensor]. Using lim e _>.o l/{x + ie) = V/x — iit5(x) (with 
V principal value), the thermal CP potential of an isotropic ground-state 
atom is given by 



77(r) =^k B T J2 'd«(4)TrG (1) (r,r,4) 

|d 0fc | 2 ^Ar( W/cj T B )TrRcG( 1 )(r,r,^). 



m=0 
k 



(7) 



The corresponding force F(r) = — VC/(r) agrees exactly with the BS result 
for the force at initial time on a ground-state atom in the perturbative limit, 
cf. Eq. (25) of Rcf. 8. Cf. this reference for details of the derivation and the 
dynamical and nonperturbative generalisation of this result for arbitrary 
(incoherent) initial-state preparation of the atom. The first term of Eq. (7) 
is the non-resonant force, while the second term is a resonant term due to 
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absorption processes, which are pure non-equilibrium effects. Note that such 
processes will unavoidably heat the atom, so that our result is only valid at 
initial time when the atom is still in its ground state. When the atom has 
become thermalized, the resonant term vanishes and the expression reduces 
to a pure Matsubara-type sum similar to the non-resonant expression. 

In the APS case, 5 the body is assumed to be held at a uniform tem- 
perature Ts different from that of the environment, Te, such that the total 
system is in a stationary non-equilibrium. The atom is assumed to be at 
zero temperature in the sense that the thermal energies fceTg, A;bTe are 
much smaller than the energies ftwk necessary to excite the atom. This 
means that terms N(w, TE)Ima(w), which are proportional to the occu- 
pation numbers of photons at the atomic frequencies, cf. Eq. (6), can be 
neglected. The density matrix p obeys the kinetic-like equation emerging 
in the Keldysh method: 15 

p(r,t;r',t') = p {r,t,r',t') + ^- j d 3 r x J dt x J ' d 3 r 2 J dt 2 [p (r, t; n, h) 

■ n j4 (ri,ti;r 2 ,t 2 ) ■ G j4 (r 2 , t 2 \ r', t') 
+ Go(r,i;ri,ti) • n R (ri,t 1 ;r 2 ,t 2 ) ■ p(r 2 , t 2 ; r', t') 

- G (r,i;ri,*i) ■ ni 2 (ri,ii;r2,*2) • GA(r2,*2;r',i')/(A«ofi)]; (8) 
G(r, t; r', t') = G (r, t; r', t') + ^ Jd s n j dt x Jd 3 r 2 J dt 2 G (r, t; n, h) 

■n fl (r 1 ,ti;r 2 ,t 2 )-G(r 2 ,t 2 ;r / ,t')- (9) 

Here, the index corresponds to free fields; Ga = G^ is the 
advanced Green tensor; and II are polarisation operators defined 
in terms of fluctuating currents IIi 2 (r, i; r', t') = —i(j(r',t')j(r,t)) T , 
n 2 i(r,t;r',0 = -i{j(r, i)j(r', t')), and U R (r, t; r', t') = [U 12 (r, t; r', t') - 
n 2 i(r, t; r', t')]@(t — t'), TIa = n^, where & is the Heaviside function. If 
the body was in equilibrium with the environment, the solution of Eq. (8) 
would be given by Eq. (3). To account for effects due to the temperature 
difference between the body and the environment, we now solve Eq. (8) it- 
eratively starting with the equilibrium solution (3). The first approximation 
reads: 

p( 1 )(r,r',w) = -2^Vo^ ^ (^i ^ ,w)ImG(r,r', W )- M 2 J d 3 ri d 3 r 2 G (r,n,cj) 

• {N(T E ,u)U2i(T U T2,u) - [N(T E ,u) + l]n 12 (r 1 ,r 2 ,w)} • G*(r 2 ,r». 

(10) 

Here, we have used Eq. (9), the free-field density matrix p (r,r',td) — 
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— 2iH/j,oN(u>, TE)ImGo(r, r', u>) and the relation 15 11^ — 11^ = II12 — II21 
for the polarization operators. Resubstituting Eq. (10) into the r.h.s. of 
Eq. (8), we obtain a result of the form (10) with Go i-> Go + Go -Ilfl- Go- 
Repeating the iterations many times and performing the sum via a Born 
series G = Go + Go-n^-Go + Go-rifl-Go-IlR-Go + . . ., we arrive at 



P(r,r',w) = -hfi 



2iN(T E ,u)ImG(r,r',u)) + J d 3 nd 3 r 2 G(r, ri,w) 



■{JV(T B ,w)n 2 i(ri,r 2) w) - [N(T E ,u) + l]n 12 (r 1; r 2 , w)} • G*(r 2 , r»] . 

(11) 

As stated, the body is assumed to be at local thermal equilibrium with 
temperature T$. This allows us to implement the fluctuation-dissipation 
theorem to calculate the polarisation operators: 



n 12 (w, r, r ') = _ !^f£ u 2 N(u, T s )Ime(r, cu)S(r - 



n 2 i(w,r,r') = - 



2tt 

^ u 2 [N(u, T s ) + l]Ime(r, u)6(r - r')- 

Z7T 



(12a) 
(12b) 



Inserting Eqs. (12) and (11) into Eq. (1), we find: 
' 2ir 



U(v) = 



dujuj 2 [2N{uj,T E ) + l]ImTr a(u) ■ G (1) (r,r,w) 

+ Rc 7TT2 / duj ^2 I rf 3 sIme(s,^){iV( W ,TB)[7V(w,rs) + l] 
Jo c J 

~[N(uj,T e ) + l]N(oj,T s )}Tr[a(uj) ■ G(r,s,w) • G*(s,r,w)] . (13) 

The first term of Eq. (13) describes the equilibrium CP force, cf. Eq. (4) 
above, the second term corresponds to the case when the medium is not 
in equilibrium with the field. If the medium is a homogeneous body of 
permittivity e(uj) occupying a volume Vs, then the latter leads to just the 
APS result for the non-equilibrium force F neq (r) = — V(7 neq (r) [see Eqs. (7) 
and (9) of Ref. 5]: 



him f°° ^ f 



d d s 



1 



e huj/k B T s _ I e hui/k B T E _ I 

x Ime(w) { V'ReTr [a(w) • G(r,s,w) • G*(s, r', w)]} r , =r 



(14) 



3. Thermal CP potential on a particle in uniform 
temperature environment 

To illustrate the non-equilibrium effects, we apply the general theory of 
the previous section to specific scenarios. We will concentrate on a non- 
equilibrium between the particle and the electromagnetic field and consider 
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the potential (7) of a ground-state particle in an environment of uniform 
temperature. 



3.1. Planar systems 

The trace of the Green tensor at a distance z to the right of a half-space 
with reflection coefficients r s ,r p for s,p polarization is 16 



Tr/ s i f°° qdq 



2/3 2 



£ 5 

(T—S,p 



(15) 



with /3 = \J us 2 / c 2 — q 2 . The integral over transverse momentum q naturally 
separates into a propagating part q < uj / c and an evanescent part q > lu/c. 

We use Eq. (15) to calculate the force on a ground state LiH molecule 
outside a gold half-space at T — 300K. The result is striking (Fig. la): the 
evanescent part almost exactly cancels the non- resonant part, and the prop- 
agating part is spatially oscillating and dominates in the retarded regime. 




Fig. 1. (a) Components of the thermal CP force F z = —dU/(dz) on a ground state LiH 
molecule outside a gold half space at 300K. (b) Enhanced potential from cavity of width 
a = Afco (solid line). The dashed line indicates the potential outside a single half-space. 



Unfortunately, the spatially oscillating propagating force component is 
very weak outside a half-space. We have investigated a scheme to enhance 
the amplitude of the oscillating potential by fine-tuning the width of a 
planar cavity to exactly one wavelength of light which resonates with the 
dominating molecular transition. 12 Figure lb shows that the scheme works 
in principle, but the enhancement factor thus achieved is not enough to 
bring the oscillations into a regime which is likely to be observable. The 
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reason for this is primarily that the enhancement factor scales with the 
logarithm of the Q-factor of the cavity (shown analytically and numerically 
in Ref. 12) which strongly limits the potentiality of such a scheme. 



3.2. Cylindrical cavity 

Another candidate geometry is a molecule situated inside a cylindrical cav- 
ity of radius R. At certain specific radii, resonances like that for the planar 
cavity are expected to occur, where the two-dimensional mode confinement 
holds the promise that resonant force components may be much stronger 
than in the planar case. As a first step towards exploiting this effect to re- 
alise strong thermal CP forces, we investigate the spatial profile of resonant 
forces inside the cylindrical cavity. Work on this problem is continuing. 

The trace of the Green tensor for points r' = r = (p, 6, z) inside a 
cylindrical vacuum cavity in an unbounded non-magnetic medium of per- 
mittivity e may, after much simplification, be written as 17 ' 18 



TrG(p, p; u)=— dt V { (r M + t 2 r N 
2lT J o 



r j^) + j; 2 (<^) 



r N ^J^x)}. (Hi) 



Here, x — g\J\ — t 2 , g = kR and k = lj/c. The dimensionless ra- 
dial co-ordinate is <j> = p/R and the integration variable t is the di- 
mensionless momentum component along the cylinder axis, h, relative to 
fc.t The reflection coefficients read tm,n = — [Hn\x)/J n {x)]rM,N with 
r M ,N = {A + B m ,n)/{A + B D ) and 

A =n 2 [x 6 - (2x 2 + g 2 )x 4 + (2g 2 + x\)x\x 2 - g 2 x% (17a) 

Bm =g 2 x\x 2 [eh\x 2 - (hij 2 + eh\h 2 )xix + h 2 j 2 xl], (17b) 

B N =g 2 x 2 1 x 2 [£h 2 1 x 2 - {eh\j 2 + h 1 h 2 )x 1 x + h 2 j 2 x\], (17c) 

B D =g 2 x\x 2 \eh\x 2 - (e + \)h{j 2 x\x + j 2 xf\ (17d) 

[xi = gVs - t 2 ; x 2 = x; jj = J' n {xj)/ J n {xj); hj = H^'(xj)/H^\xj)]. 
In the perfectly conducting limit |e| — > oo (at nonzero lj) wc find 

H { n ] '(x) H£\x) , 1C , 

rM ^-^r ; rw ^w (18) 



^For the non-resonant term the dimensionless substitution variable t = h/k cannot be 
used since k = for the zeroth Matsubara term. We use i = hR = gt in this case. 
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For the perfectly conducting cylinder, thus, the resonant radii for radiation 
of a given frequency are R n i{u) and R' ni (uj) with R^j (w) = cj^j/uj wherein 
j ni and j' ni are the ith zero of J n (x) and J' n (x) , respectively. When e < oo 
the resonances move away from these values. Further details and analysis 
will be reported elsewhere. 




700 



Distance from axis (jim) 



Fig. 2. The thermal CP potential on a LiH molecule in an infinitely thick cylindrical 
gold cavity at T = 300K. The radius is R = 1.5Ru « 618(im. 

We plot the potential for the example of R = 1.5J?n for ground-state 
LiH in Fig. 2. For technical reasons, in the cylindrical geometry splitting 
the resonant potential into propagating and evanescent parts is no longer 
natural and straightforward. As a benchmark we have checked that the 
potential outside a plane is regained close to the cylinder boundary in the 
limit of large cylinder radius R ^> 2ir/kiQ. 

Figure 2 shows clearly that the peculiar traits observed for the resonant 
potential outside a half-space, depicted in Fig. la, are present also in the 
cylindrical cavity as one would expect. The resonant potential once again 
almost cancels the non-resonant term close to the surface giving a resulting 
attractive force in the near zone which is dramatically reduced compared to 
the non-resonant term alone. As before the retarded regime is dominated 
by oscillating behaviour. However, we also observe qualitative differences 
from the planar case, e.g., the potential minimum closest to the wall is no 
longer the deepest one as is the case in a planar cavity. 12 

4. Summary 

We have demonstrated that the recent complementary theories for the ther- 
mal CP force between a ground-state particle and a body at non-equilibrium 
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between either the particle or the body and the electromagnetic field may 
both be obtained using the Keldysh formalism. Applying the results to 
planar and cylindrical geometries of uniform temperature, we have found 
that even a ground-state particle is subject to resonant force components 
which in both geometries strongly cancel the well-known resonant force in 
the nonretarded regime and lead to spatially oscillating forces for retarded 
distances. The latter are expected to be much enhanced for the cylindri- 
cal cavity, opening the perspective towards molecule guiding with thermal 
photons. 
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